Electron-phonon interaction enhanced by antiferromagnetic and superconducting 
fluctuations in cuprate oxide superconductors 
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An electron-phonon interaction arising from the modulation of the superexchange interaction by 
phonons is studied within the theoretical framework of Kondo lattices. It is relevant in strongly corre- 
lated electron liquids in cuprate oxide superconductors, which lie in the vicinity of the Mott-Hubbard 
metal-insulator transition. It is enhanced by antiferromagnetic and superconducting fluctuations, 
which are developed mainly because of the superexchange interaction. When the enhancement of the 
electron-phonon interaction is large enough, it can explain the softening of phonons and kinks in the 
quasiparticle dispersion in cuprate oxide superconductors. However, the superexchange interaction 
itself must be mainly responsible for the formation of Cooper pairs. 

PACS numbers: 74.20.-z, 71.38.-k, 75.30.Et 



I. INTRODUCTION 

It is an important issue to elucidate the mechanism of 
high-T c superconductivity occurring in cuprate oxides. 1 
It was shown in previous papers^ published in 1987 that 
the condensation of c?7-wave Cooper pair bound by the 
superexchange interaction can explain observed T c ; T c for 
d'y wave are definitely much higher than those of other 
waves, as long as the on-site repulsion is so strong that 
cuprate oxides with no dopings might be Mott-Hubbard 
insulators. On the other hand, two observations, the 
softening of phonons^ 6 *^ and kinks in the quasipar- 
ticle dispersionj2ii2iAiii2iiiii4 imply the relevance of an 
electron-phonon interaction. In particular, the softening 
of phonons is evidence that an electron-phonon interac- 
tion is strong in cuprate oxides. One may argue that it 
must be responsible for high-T c superconductivity or, at 
least, it must play some role in the realization of high-T c 
superconductivityiiSiii 

Various experiment o 15 ' 16 i 17 ' 18 i 19 ' 20 i 21 ' 22 imply or show 
the opening of anisotropic pseudogaps at temperatures 
above T c in quasiparticle spectra of the so called optimal 
and under-doped cuprate oxide superconductors. One 
may argue that kinks are caused by or are closely re- 
lated with what cause pseudogaps, rather than phonons. 
Antiferromagnetic (AF) and superconducting (SC) fluc- 
tuations are developed in cuprate oxide superconduc- 
tors; not only antiferromagnetism and superconductiv- 
ity themselves but also the development of their fluctu- 
ations are mainly caused by the superexchange interac- 
tion. It was shown in a previous paper 2 - that large life- 
time widths, which are mainly caused by well developed 
d7-wave SC fluctuations so that they are anisotropic or 
wave-number dependent, are responsible for anisotropic 
pseudogaps; spectral weights around the chemical po- 
tential are swept away because of large life-time widths. 
Small kinks appear in calculated spectra^ but they are 
too small to explain observed kinks. It is difficult to ex- 
plain kinks by AF or SC fluctuations. It is reasonable 
that not only the softening of phonons but also kinks are 
caused by an electron-phonon interaction. 



Doped holes mainly go into O ions. This implies that 
the local charge susceptibility of 3d electrons on Cu ions 
is much smaller than that of 2p electrons on O ions and 
charge fluctuation of 3d electrons can never be devel- 
oped. It is quite unlikely that the conventional electron- 
phonon interaction, which directly couples with charge 
fluctuations, plays a crucial role in cuprate oxide super- 
conductors. 

A necessary condition for a relevant electron-phonon 
interaction is that it can work even in strongly corre- 
lated electron liquids in the vicinity of the Mott-Hubbard 
metal-insulator transition or even when charge fluctu- 
ations are significantly suppressed. One of the most 
plausible ones is an electron-phonon interaction arising 
from the modulation of the superexchange interaction 
by phonons. It was pointed out that it plays a role 
in phonon-assisted multi-magnon optical absorption— It 
can couple directly with AF and SC fluctuations, so that 
it can be substantially enhanced by AF and SC fluctua- 
tions. One of the purposes of this paper is to show that 
the electron-phonon interaction arising from the modula- 
tion of the superexchange interaction by phonons is rel- 
evant, at least, in cuprate oxide superconductors where 
AF and SC fluctuations are substantially developed. 



II. FORMULATION 
A. Electron-phonon interaction 

It was shown in a previous paper— that Gutzwiller's 
quasiparticle band 26 lies between the lower and upper 
Hubbard bands^l in metallic phases in the vicinity of the 
Mott-Hubbard transition. Gutzwiller's quasiparticles are 
responsible for metallic properties. The superexchange 
interaction arises from the virtual exchange of pair exci- 
tations of electrons across the lower and upper Hubbard 
bands. As long as the Hubbard splitting is significant, 
therefore, it works between Gutzwiller's quasiparticles. 

When we follow previous papers22i22^S and we ignore 
nonzero bandwidths of the lower and upper Hubbard 
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bands, it is straightforward to show that the virtual ex- 
change process gives the following exchange constant be- 
tween nearest-neighbor ith and jth Cu ions: 
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where the 3d levels of the ith and jth Cu ions are de- 
noted by tdi and e d j, the 2p level of the [ij]th O ion 
between the two Cu ions by e p u,j] , and the hybridization 
energies between the Cu ions and the O ion by V^m and 



Vj t [ij], respectively. When we put e dl — > e d , e dj - 
^'tfa'] = = V , we obtain a wellknown one: 
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The variation of Jy is given by 
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When we take the x- and j/-axes along Cu-O-Cu bonds, 
variations of e d i, e p [yi and V^ryi are given by 

Ae e2 i = J 4 d [e a: -(ui, x+ -Ui j:(; _)+ej / -(ui,j /+ -Uj, !/ _)] , (2.4) 
Ae p[lj] = A p [ eij ■ (uj - Uj)] , (2.5) 



AVij^] + AVj-jjj] = Ay [e^ • (u< - Uj)] 



(2.6) 



to linear order in displacement of ions, with A d , A p and 
Ay being constants, U; the displacement of the ith Cu 
ion, Ui ^ s that of an O ion on the adjacent s = + or 
s = — side along the £-axis of the ith Cu ion, e x — (1, 0), 
e y = (0, 1), and ey = (R; — Rj)/|Rj — R,-|, with Rj the 
position of the ith Cu ion. 

Displacements of the ith Cu and the [ij] th O ions are 
given by 



Hv, 
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with Rryi = (l/2)(Ri +Rj), M d the mass of Cu ions, 
M p the mass of O ions, 6>q and & A _ q annihilation and 



creation operators of phonons with polarization A and 
wave vector q, u^q energies of phonons, e Aq unit polar- 
ization vectors, and N the number of unit cells. The q 
dependence of v d ,\q and v p ,Aq can play a crucial role. For 
example, Vd,\q = and i> p ,Aq = 0(1) for breathing modes 
that bring no changes in adjacent Cu-Cu distances. 

The electronic part can be well described by the t-J 
model on a square lattice: 

H = e d ^24 a .d ir7 - 1 ^2t ij 4 a d j( r 

i<j ija 

(ij) 1 

with the summation over (ij) restricted to nearest neigh- 
bors, 



(2.10) 
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with a x , Oy and a z the Pauli matrices, and ni a = d\ a di a - 
An infinitely large on-site repulsion, £/oo/l^(y)l ~~ * 
is introduced to exclude any doubly occupied sites. 

According to Eq. I|2.4[) . there are two types of electron- 
phonon interactions. Define an operator by 

?Mq) = 5$>r(q')[S (q'+§q)- S (-q'+|q)] , (2.11) 



q' 



(2.3) with 



s (q) = -Ts E l^cM-w-Vw/f ( 2 - 12 ) 



with cr = (<7 X , (jy,cr z ). They are given by 



iC p £ 



p,Aq 



&I-q + ^q 



^ yj / 2NM p (jJ\q 

U<l) E ^(|q)Pr(q), 



(2.13) 



r=s,d 

^ ^tq + ^Aq 



y/2NM d U X q 



E *?r(q)Pr(q), 



r=s,d 



with 



8A d V 4 



(e d + U - e p ) 3 



4A p V 4 



e d + U-e p U 



(e d + U-e p y 



e d + U-e p U 



2A V J 



q x . fq x a\ q y . (q y a\ 
— sin — — + — sin -2— 
a \ 2 J q \ 2 J 



(2.14) 



(2.15) 



, (2.16) 



(2.17) 



3 



% (q}=2 -sm(— j sm( — j , (2.18) 



with 



?7 s (k) = cos(fc K a) + cos(fc y a), (2-19) 



T]d(k) — cos(k x a) — cos(fc y a), (2.20) 

with a the lattice constant. Here, we consider only lon- 
gitudinal phonons or e^ q = {q_x,<lyi<lz) / 'q is assumed. 



B. Theory of Kondo lattices 

We follow the previous paper^- to treat the infinitely 
large Uoc, where a theory of Kondo lattice is developed. 
A renormalized single-site approximation (SSA), which 
includes not only all the single-site terms but also the 
Fock term AE(k) due to the superexchange interaction, 
is reduced to solving the Anderson model with the in- 
finitely large on-site repulsion U<x,. The self-energy of 
the Anderson model is expanded as 



S CT (i£„) = S(0) + (1 - 4> 7 )is r , 



■5^(l-^)A/i 



, (2.21) 



with A/v infinitesimally small spin-dependent chemical 
potential shifts. Note that </> 7 = (j> aa . The Wilson ratio 
is defined by W s = s /</> 7 , with <f> s = 4> aa - 4> a - a . For 
almost half filling, charge fluctuations are suppressed so 
that 4> c — 4> aa + 4> a -a *C 1. For such filling, 7 ^> 1 



so that 



quasiparticles is given by 



For such filling, ^ 7 
2<f> 1 or W s — 2. The dispersion relation of 



07 



■E(0)+AE(k)-Ai , 
(2.22) 



with /i the chemical potential. 

The spin susceptibility is given by 



Xs(«^,q) 



27r s (iwi,q) 



1- [^J(q) + Uoo] TT s (iu}i,q)'' 



(2.23) 



with 7r s {iuji , q) the irreducible polarization function in 
spin channels and 



J(q) = 2j7? s (q). 



(2.24) 



The function -ir s (iuji, q) is divided into single-site Tr s (iuji) 
and multi-site Air s (icui , q) in such a way that ir s {iuui , q) = 
tt s {iuji) + ATT s (iuji, q). In Kondo lattices, local spin fluc- 
tuations at different sites interact with each other by an 
exchange interaction. Following this physical picture, we 
define an exchange interaction I s {iwi , q) by 



Xs(«w;,q) 



(2.25) 



2jT s (iu)i) 



1 - UooTTs^Wl) 
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the susceptibility for the mapped Anderson model. Then, 
we obtain 



I s (iwi,q) = J(q) + 2J7^A7r s (iw/,q). 



(2.27) 



The main part of 2U^ c Air s (iuji, q) is an exchange interac- 
tion arising from the virtual exchange of pair excitations 
of quasiparticles. 

When the Ward relational is made use of, the ir- 
reducible single-site three-point vertex function in spin 
channels, X s (ie n , ie n +iuji; iu>i), is given by 

UooXsiienJSn + iajfjuji) = 2<p s /x s (iuji), (2.28) 

for |e„| — > +0 and \uii\ — > +0. We approximately use 
Eq. fT2"8l for |e„| < k B T K and \u x \ < k B T K , with T K 
the Kondo temperature defined by 



k B T K = [1/X S (0)] 



T-+0 ' 



(2.29) 



The so called spin-fluctuation mediated interaction, 
whose single-site term should be subtracted because it 
is considered in SSA, is given by 



1 

4 

with 



and 



■1.k!\ ji^,) 5 F(»w,,q) = ^i/;(^,q), (2.30) 
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-I*{iui,q) = - 



Tl s (iu} h q)x s (iuJi) 
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(2.32) 



Because of these equations, we call / s (iw;,q) a bare ex- 
change interaction, I*(iuji, q) an enhanced one, and <p s an 
effective three-point vertex function in spin channels. In- 
tersite effects can be perturbatively considered in terms 
of F(ioji,q), I s (iuji,q) or I*(ioji,q) depending on each 
situation. 

The enhanced one is expanded as 

r s (iu h q) = J * + 27 1 *ry s (q) + 27 2 *ry s2 (q) + • • • , (2.33) 

with 

Vs2(q) = cos[(k x +k y )a] + cos [(k x -k y ) a] . (2.34) 

The nearest-neighbor I£ is mainly responsible for the 
development of SC and charge bond-order (CBO) 
fluctuations^ Because contributions from < k B Tx 
are the most effective, we ignore its energy dependence. 
An effective SC susceptibility is calculated in the ladder 
approximation: 



r (SC) 



XrJdWuuq) = 



(2.35) 
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for r = d wave or dry wave, with 



r (SC) 



(iw/,q) 



nk 
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ie„-£(k+±q) 



-ie n - iui - £(-k + iq) 



(2.36) 



Only c?7-wave SC fluctuation are considered in this paper 
because T c of d'y wave are definitely much higher than 
T c of other waves An effective CBO susceptibility is 
similarly given by 



(CBO)/- \ 



'(^,q) 



i + f/r^4 CBO) (^,q)' 

for r = s, p and c?7 waves, with 



(2.37) 



'(iwj,q) 



ie„-£(k-|q) 



fe„ + icj; - £(k + iq) 



The form factors of p waves are defined by 
77a; (k) = v^sin^o), VyO*-) = V2sin(k y a). 



(2.38) 



(2.39) 



According to Eq. H2.35(l . T c of ^7 wave superconduc- 
tivity are given by 



i + |/r^>< sc, (o,o) = o. 



It was shown in the previous papers^ that 
f \I*\W^ ~ 100 meV 



(2.40) 



(2.41) 



is needed in order to explain observed T c . The Wilson 
ratio is as large as W s ~ 2 in SSA, and the superexchange 
interaction is as strong as J = —(100-150) meV in ac- 
tual cuprate oxide superconductors. Then, it follows that 
|l J rl^ s 2 Z 400 meV; 1^*1 > \ J\. The theory published in 
1987 has a drawback that it gives too high theoretical T c 
to explain observed T c m£ The mass enhancement factor 



and the effective three-point vertex function are renor- 
malized by AF, SC and CBO fluctuations, so that W s 
that is the effective three-point vertex function divided 
by the mass enhancement factor is also renormalized by 
the fluctuations. We argued in the previous paper— that 
the renormalization of W s is substantial so that W s ~ 1 
or W s < 1; a phenomenological argument also implies 
that W s = 0.7-1 had better been used in order to explain 
quantitatively T c and T-linear resistivities. For example, 
Eq. (|2.41() can only be satisfied for such small W s . We 
follow this argument. Taking the renormalization into ac- 
count, we regard W s as a phenomenological parameter; 
we assume 



W„ = 1 



(2.42) 



in this paper. 



C. Renormalization of phonons 

The Green function for phonons is given by 

2w Aq 



D x (iuJi,q) 
with 



(iivi ) 2 - + 2uj Xci Aujx {iuji , q) ' 



Acj A («wz,q) 



2M p UJ\q 



(2.43) 



(2.44) 



Because phonons are renormalized by AF, SC and CBO 
fluctuations as well as charge density fluctuations, we 
consider four processes shown in Fig. ^ S(iLUi,q) — 
S s (ioJi,q) + Ssc{ iuJi,q) + 5 CBO (^/,q) + S c (ioj t ,q). The 
expression <|2.11|) for Vr (q) is convenient in treating cou- 
plings with AF fluctuations. Other expressions, which 
are convenient in treating couplings with SC and CBO 
fluctuations, are shown in Appendix. When Eqs. ijAlfl . 
(|A3(I . and l|A6[) are made use of and only the parts of 
r = s in Eqs. (|2.13|l and (|2.14|l are considered, it follows 
that 



J 



S s (iui,q) = ^y A 2 (q)^^^77^(q')x s (^/ + ^',q'+ 2 -q)x s (^* CJ i''^q' + 5q) > 

;'q' 



(2.45) 



5 S o(?wz,q) 



^X 2 (q)^E 
I'd' 



Xf (iuJi + uj v , q' + |q) xf (-«"!' , - q' + |q) 



-n s d ° (iui + iu) V , q' + ±q) nf (-iu) V , -q' + ^q) 



(2.46) 
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FIG. 1: Four processes of the renormalization of phonons. A solid line stands for an electron, a broken line for a phonon, a 
wavy line for the superexchange interaction J, and a solid circle for the effective vertex function <j> a - Hatched parts in Figs.^a), 
(b) and (c) stand for AF (spin), SC and CBO fluctuations, respectively, and a hatched part, including an internal electron line, 
in Figs. rTTd) for the vertex function Z(ie n , iuii; k, q). The contribution from Figs. Q^d) is larger than those of Figs.^a), (b) 
and (c); two fluctuation-lines with different q or two susceptibilities of different q appear in the convolution form in Figs.^a), 
(b) and (c) so that their contributions can be large only when fluctuations are developed in a wide region of the momentum 
space. 



N 

r Z'q' 



Xr BO (iwj + iu v , q' + |q) Xr B ° {-m> , -q' + |q) 
ttT° (t^i + iwi/ , q' + |q) t^T° , -q' + H 



(2.47) 



S c (iu>i, q) 



4 2 



1 



nker 



ie n -£(k) ie n +iui- £(k+q) : 



with 



F A (q)=7? s (q) 



CpVp^qTjs (|q) + C d Vd,XqJ M p /M d 



(2.48) 



(2.49) 



In Eq. H2.48fl . Z(ie n) icj;;k, q) is the vertex function in the charge channel. It is also enhanced by AF, SC and CBO 
fluctuations; Z(ie n , iwr, k, q) = Z s (ie n , iu>i; k, q) + Z sc (ie„, k, q) + Z CBO (ie„, iLu t ; k, q) H , with 



Z s (i£ n ,iui l ;'k, q) 



N 



AT S q' + |q) AT S (-io;;/ + -q' + |q) 



j'q' 



ie n + iw'j - £(k + q' + ±q) 



(2.50) 



^sc(jen,iwz;k, q) 



1 k B T \ 

2 AT ^ 

J'q' 



% (k 



|q) Vd (k - iq' + |q) 



Kf (i^.q'-Iq)^, 



«w ; ,q 



+ icJ; - £(— k + q' 



(2.51) 



with 



^cBo(iffn 5 ^;k,q) 



l fciiT 
"2 AT 



EE^C'+^ + Sl)^ (k+|q' + |q) 



J'q' r 



^ BO (io;,,,q' + ±q) A^ BO (-fag + z^, q' + |q) 
ie n + iuj[ - £(k + q' + |q) 

I 



(2.52) 



AT s (iw;,q) = 1 ! 



1 



l- 3 /(itj;,q)x s («Wi)' 



l + fjfW?*§°(iw,,q) 



(2.53) 
(2.54) 



^ BO (^,q) 



l + lJTWj^fiwi.q) 



(2.55) 



Here, Eqs. (531, and (5U are also made use of; 

zero-th order terms in /* are subtracted in Eqs. I|2.46|l . 
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(|2.47|) , I2.54|l and (|2.55|) to avoid any double counting. 



When we put 



III. APPLICATION TO CUPRATE OXIDE 
SUPERCONDUCTORS 

A. Softening of phonons 

Because it is not a purpose of this paper to study 
phonon modes themselves, we study only issues that can 
be clarified without calculating them. First, we consider 
the softening of the so called half breathing mode of O 
ions with Qx = (±7r/a,0) or (0,±7r/a). Because 



0, 



1. 



(3.1) 



we have to consider only the electron-phonon interaction 
given by Eq. 12.1311 . 

Because doped holes mainly go into O ions, one may 
argue that 3c? levels of Cu ions are deeper 2p levels of 
O ions, e<2 < e p . However, this argument disagrees with 
what band calculations^iMiS predict, — e p — 1 eV. 
The observation does not necessarily mean that < e p , 
but it simply means that the local charge susceptibility of 
3d electrons is much smaller than that of 2p electrons, as 
is discussed in Introduction. Because it is unlikely that 
band calculations give such a bad prediction on relative 
positions between the 3d and 2p levels, we follow what 
band calculations predict. When we use 



V=1.6eV, e d 



e p = 1 eV, 



U = 5 eV, 



(3.2) 



following the previous paper^S Eq. 12.21) gives J — 
—0.27 eV. This is about twice as large as the experi- 
mental one of J = —(0.10-0.15) eV. This discrepancy 
is resolved when nonzero bandwidths of the lower and 
upper Hubbard band are considered^ Because of this 
reduction, we assume a half of C p given by Ea. (|2.15|l : 
C p = 0.1 x Ad instead of C p = 0.22 x Ad- When we take 
A d ~ 5 eV/A, it follows that 



C p ~ 0.5 eV/A, Y x (Qx) ^ 1 eV/A. 



(3.3) 



Two susceptibilities appear in the convolution form in 
Eqs. (E3SJ), (|2~4^|) and fjtgfi. Unless AF, SC and CBO 

fluctuations are developed in a wide region of the mo- 
mentum space, the convolutions cannot be large. We 
assume that S c (iu>i, q) given by Eq. (|2.48l) is dominant in 
S(iuJi,q). When only S c (iuji,q) is considered, 



3\ T , 



S^,q)^-V^ 4 (Z 2 ) 



x) 



2k B T K 



(3.4) 



with (Z 2 ) an average of Z 2 (ie n , iwi;k, q). When we as- 
sume W s — 1, the softening at X point is given by 

Acj A (^A qx >Qx) -0.01 (Z 2 ) (1 ° 3 meV. (3.5) 



WAQx — k B T K ~ 10 2 meV, 



it follows that 



A WA ( WAqx ,Q x ) ~ (Z 2 ) meV. 



(3.6) 



(3.7) 



When AF, SC and CBO fluctuations are not developed, 
{Z 2 ) < 1 and the softening must be very small. When 
they are well developed so that (Z 2 ) might be as large 
as (-Z 2 ) ~ 10, we can explain the observed softening as 
large as 



Aw A (wAq x ,Qx) - 10 meV. 



(3.8) 



When other contributions are considered in addition to 
S c (iu)i,q), (Z 2 } can be smaller than 10 to explain the 
observed softening. 

No softening occurs for q = because fj s (0) = 0. 
When q goes from V point to X point, the softening 
must increase first but it is unlikely that the softening 
is the largest at X point. Because v dt \Q x = 0, the 
electron-phonon interaction described by Eq. I|2.14|) van- 
ishes. This implies that the softening cannot be the 
largest at X point along T-X line. In actual, several 
experimental data imply that the softening is the largest 
for q a little different from Qx along T-X lineA^ 

No softening cannot occur either at M point or for 
the breathing mode of O ions with Qm = (±7r/a, ±7r/a) 
because v d ,\Q M = and J? s (|Qm) = 0. It is interesting 
to confirm this prediction. 



B. Kinks in the quasiparticle dispersion 

A process corresponding to Fig. 1(d) renormalizcs 
quasiparticles. The self-energy correction is given by 

J-AE(ie n ,k) = --fp Y] gl(ie n , iur, k, q)D x (iuj[, q) 

: : . (3.9) 



is n + iuji - £(k + q) ' 



with 



g\{ie n , iun; k,q) = 



v /2AfpWAq 

xZ(i£„,iu;;;k, q). (3.10) 



It is likely that the contribution of Fig. Cfd) dominate 
those of the other three, Figs.^a)-(c). In such a case, 



g\{ie n , iwjjk, Qx) ^ \J 2k B T K \&u\{uj\c L , Qx)\- (3.11) 

Here, Eas. l|2.44|) and l|3.4|) are made use of. When the 
experimental value (|3.8|l is used, we obtain 



g\(ie n ,iuji; k, Qx) ^ 45 meV. 



(3.12) 
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This is large enough for optical phonons to cause kinks 
in the quasiparticle dispersion. 

Two types of kinks are observed^ The renormaliza- 
tion by phonons can explain one type of kinks observed in 
both normal and SC phases. However, it is difficult to ex- 
plain the other type of kinks observed only in SC phases; 
low-energy AF and SC fluctuations are suppressed when 
SC gaps open. 

C. Cooper-pair interaction 

The phonon-mediated pairing interaction is given by 

V ph (q; k) = -2^(0, 0; k, q)/^ Aq . (3.13) 

Its average over k on the Fermi surface is expanded in 
such a way that 

(V ph {q;k)) =V a + 2V 1 7 ls {q) + 2V 2 r )s2 (q) + • • • . (3.14) 

No softening at T and M points implies (V^,^(0;k)) = 
and (V^(QM;k)) = o, so that 

V + 2Vi + 2V 2 ~ 0, Vq - 2Vx + 2V 2 ~ 0. (3.15) 

Because of Eq. i[3~T2|l . (V ph {Q x ; k)) ~ -40 meV or 

V - 2V 2 ~ -40 meV. (3.16) 

Then, we obtain 

V ~ -20 meV, Vi ~ meV, V 2 ^ 10 meV. (3.17) 

The interaction Vi between nearest neighbors should 
be included in addition to |7* W$ in the theory of thf- 
wave high-T c superconductivity. When Eq. (|2.40() is ex- 
tended to include Vi, T c are determined by 

1 + (Vriy s 2 + Vi) 7r^ sc, (0,0) = 0. (3.18) 

Although the q dependence of (V^(q;k)) in the whole 
Brillouin zone is necessary to estimate Vi accurately, 
we can conclude that |Vi| must be much smaller than 
ll/j^W^ 2 ~ 100 meV, which is needed in order to explain 
observed T c m£. 

There are various branches of phonon modes beside the 
mode discussed above. The virtual exchange of phonons 
that do not become soft cannot give a significant pairing 
interaction. 

IV. DISCUSSION 

Following Barnes?2& we can map the t-J model to the 
so called auxiliary-particle t-J model: Ht-j = T'~ 1 'H'P, 
with 

^ = 6 E ( e * e * + c W + 4i°ii - X ) + e d 4° Cl ° 

i ia 
ija (ij) 



with S being an arbitrary constant^ and 

S, = ^W^^)i^- (4-2) 

a/3 

Two kinds of auxiliary particles, which correspond to 
empty and occupied sites in the original t-J model, are 
introduced; e] and c| CT are their creation operators. We 
call them e and c particles in this paper. The projection 
operator V restricts the Hilbert space within 

Qi = eje; + cj T c iT + c- x c a = 1 (4.3) 

for any i site; no empty or multiply occupied sites are al- 
lowed. This restriction is guaranteed by the conservation 
of the number of auxiliary particles such as 

[n,Qi]=0 (4.4) 

for any i, or local gauge symmetry. This symmetry is 
inherent in the auxiliary-particle model. Local gauge 
symmetry can never be broken^ Therefore, no single 
auxiliary particle can be added or removed, or no single- 
particle excitation of auxiliary particles is allowed. Aux- 
iliary particles themselves are never itinerant but are lo- 
calized; pair excitations of auxiliary particles are itiner- 
ant. Fermionic pair excitations of auxiliary particles cor- 
respond to electrons in t-J model. Two ways of statistics 
are possible: fermionic e and bosonic c particles, and 
fermionic c and bosonic e particles. The model with 
bosonic e and fermionic c particles is often called the 
slave-boson t-J models The mean-field (MF) theory for 
H, instead of V~ lr H.T ', and its more or less improved 
theories, which include gauge fluctuations, assume the 
breaking of local gauge symmetry, and they treat single- 
particle excitations of itinerant auxiliary particles. Such 
theories are never relevant to study dynamics of electrons 
in the original t-J model; condensation energies derived 
in these treatments are consistent with Gutzwiller's the- 
ory and are reliable. A more precise discussion on this 
issue can be found in Appendix of Ref. |2jJ. Although 
an apparently similar theoretical development to that of 
this paper is possible when one starts from the MF ap- 
proximation for the slave-boson t-J model, it is physi- 
cally and essentially different from the theory of Kondo 
lattices; states considered in the two theories are of to- 
tally different symmetry from each other. Therefore, we 
should abstain from comparing results based on the the- 
ory of Kondo lattices with those of the MF theory of the 
slave-boson t-J models 

Note that <f> c and l/</> 7 are small parameters in the 
vicinity of the Mott-Hubbard transition. What are con- 
sidered in this paper are of leading order in both <j) c and 
l/(f> 7 , that is, order of (0 c )°(l/(/) 7 ) . 

There are two other types of electron-phonon interac- 
tions: the modulation of 3d-electron levels, ed, and that 
of the transfer integrals, Uj. The conventional one aris- 
ing from the modulation of e^, which can directly cou- 
ples with charge fluctuations, gives renormalization ef- 
fects higher order in <f> c and l/<^ 7 , so that its effects must 
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be very small. The electron-phonon interaction arising 
from the modulation of tij gives renormalization effects 
higher order in 1/4> 1 , so that their effects of the electron- 
phonon interaction must be l/0 7 times as small as those 
studied in this paper; we expect that coupling constants 
for Uj, which correspond to Ad, A p and Ay of this pa- 
per, are of the same oder of magnitude as Ay . Then, we 
ignore both of them in this paper. 

There are also pieces of experimental evidence that 
the electron-phonon interaction arising from the modu- 
lation of €d or by phonons is irrelevant in cuprates. 
Hwang, Timusk and Gu investigated life-time widths of 
quasiparticles instead of kinks; life-time widths depend 
on temperature and dopings. 1 — Their observation implies 
that kinks are large only in metallic cuprates where AF 
and SC fluctuations are well developed. The softening of 
phonons is also large in such metallic cuprates; no signif- 
icant softening is observed in over-doped cuprates— It 
is difficult to explain these observations in terms of the 
electron-phonon interaction arising from the modulation 
of e e i or t^ by phonons. On the other hand, these observa- 
tions are pieces of evidence that the electron-phonon in- 
teraction that can couple directly with AF and SC chan- 
nels is relevant and its enhancement by AF and SC fluc- 
tuations in metallic phase is crucial. 

In insulating phases, only the contributions from 
Eq. (|2.45|) remains but those from Eqs. H2.4(j[) - I|2.48[) van_ 
ish. Then, we cannot expect significant softening of 
phonons. This is also consistent with experiment— 

Various physical properties are different or asymmet- 
ric between hole-doped and electron-doped cuprates. 
Within the theoretical framework of this paper, hole- 
doped and electron-doped cuprates must be, in essence, 
similar to each other. Phcnomcnologically, AF and SC 
fluctuations are relatively more developed in hole-doped 
cuprates than they are in electron-doped cuprates. If the 
asymmetry of the fluctuations can be explained, we can 
explain that of phonon properties. It is pointed out in 
another paper— that the asymmetry of disorder can play 
a crucial role in the asymmetry between hole-doped and 
electron-doped cuprates. 

In cuprate oxide superconductors, the exchange inter- 
action arising from the virtual exchange of pair excita- 
tions of quasiparticles is less effective than the superex- 
change interaction; the pairing interaction arising from 
phonons can play no significant role. Then, the main 
pairing interaction /j* must arise from the superexchange 
interaction, which is enhanced by spin fluctuations. 

As is discussed in Introduction, the superexchange in- 
teraction arises from the virtual exchange of pair exci- 
tations of electrons across the lower and upper Hubbard 
bands. As is shown in Eqs. i|2.3U[) and (|2.32[1 . the spin- 
fluctuation mediated pairing interaction is essentially the 



same as the superexchange interaction if high-energy spin 
fluctuations, whose energies are as large as the Hubbard 
onsite repulsion U, are properly included. However, it is 
physically different from the superexchange interaction if 
only low-energy spin fluctuations are included. 

The SSA is rigorous for Landau's Fermi-liquid states in 
infinite dimensions^ so that the theory of Kondo lattices 
can be regarded as a 1/d expansion theory, with d being 
the spatial dimensionality. One may suspect that the 1 jd 
expansion theory cannot be applied to quasi-two dimen- 
sional cuprates. Any perturbative theory relies on the 
analytical continuity^ a perturbed state must be of the 
symmetry as an unperturbed state is. Normal states in 
over-doped or optimal-doped cuprates are certainly Lan- 
dau's normal Fermi liquids. Because there is no evidence 
that any symmetric change occurs between normal states 
in over-doped or optimal-doped cuprates and exotic nor- 
mal states in under-doped cuprates, we can argue that 
the analytical continuity holds so that the I jd expansion 
theory or the theory of Kondo lattices can be applied to 
exotic normal states in cuprates. 



V. CONCLUSION 

The electron-phonon interaction arising from the mod- 
ulation of the superexchange interaction by phonons is 
relevant for strongly correlated electron liquids in the 
vicinity of the Mott-Hubbard transition. It is shown with 
the help of the theory of Kondo lattices that it can be 
enhanced by spin, superconducting, and charge bond- 
order fluctuations as well as charge fluctuations. The en- 
hanced electron-phonon interaction is responsible for not 
only the softening of phonons but also kinks in the dis- 
persion relation of quasiparticles in cuprate oxide high- 
temperature superconductors. However, it can never be 
the main Cooper-pair interaction. The main one must 
be the superexchange interaction. 
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APPENDIX A: VARIOUS EXPRESSIONS FOR 
THE ELECTRON-PHONON INTERACTION 

Equation l|2.II|l is also written in another form: 
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'Pr(q) = t^t ^ ^ »7r(q') (s Q/3 - s 7(5 ) at, ,i , ,1 , »n<_ i„/i/?o! i ,,i , a (T .,i n , 

2N ^— ' ^— ' (k+iq'+iq)a ( K ^ (p-^q' + jq)! lP+ 2 q 



(Al) 



kpq' a/?7<5 



This expression is useful to obtain Eqs. (|2.45() and l|2.50[l . When we replace variables in Eq. IjAlfl in such a way that 
k+ iq' = ki + |qi, k — |q' = Pi + ^qi, and p + |q' = ki — |qi, Eq. IjAljl turns out to 



V ^ = 2^ E E ^ - Pl) ( SQ/5 ' S75 ) a [ kl + i qi + lq )ct a (p 1 + iq 1 )/3 a fp 1 -i qi + l q h a (k 1 -iq 1 )^ ( A2 ) 
kipiqi afiyS 

Matrix elements of a — S and /3 = 7 are relevant for CBO channels, and they are the same as those given by 
^f(q) = E ^(ki _ Pi)^ [ _3 Pc(ki + |qi + |q,ki - |qi)p c (pi - |qi + |q,pi + |qi) 



kipiqi 



-p s (ki + |qi + iq,kx - |qi)p s (pi - |qi + |q,Pi + ^qi)] 



with 



Pc(ki,k 2 ) = a^^k^ +aj cii a k2i , p s (k 1; k 2 ) = al^a k2l - a^ak 



2i> 



(A3) 



(A4) 



except for those given by what appear through the commutation of operators, which are not shown here. This 
expression is useful to obtain Eqs. (|2.47|) and l|2.52|) . When we replace variables in Eq. (|A2|) in such a way that 
k i + |qi = k 2 + -|q 2 , Pi + |qi = P2 + |q2, and ki - ±q x = -p 2 + ±q 2 , Eq. (jA"2|) turns out to 



Vr{ ^ = ^N E E r 'r( k 2-P2)(s Q/3 -s 7 ')a; k2+ i q2+ i q)Q «(p 2 +^ 2 )/3 a [ 



t (-k 2 + iq 2 + iq)7 (-P2 + ^q 2 )5- 



(A5) 



k 2 p 2 q 2 aflyS 

Matrix elements of 7 = —a and 5 = —(3 are relevant for singlet SC channels, and they are the same as those given by 
^f(q) = 4 E 7?r(k 2 -p 2 )i -3p|(k 2 + iq 2 + iq,-k 2 + iq 2 + iq)pi(p 2 + iq 2 ,p 2 + iq 2 ) 



k 2 p 2 q 2 



with 



+4(k 2 + ±q 2 + |q, -k 2 + ±q 2 + ±q)p 3 (P2 + |q2, -P2 + ^2) 

4(ki,k 2 ) = a klT a k2i +4 1 i4 2 T' 



pI(ki,k 2 ) = 4 lT 4a-4a4 2 T 



(A6) 



(A7) 



except for those given by what appear through the com- 
mutation of operators, which are not shown here. This 
expression is useful to obtain Eqs. I|2.46|l and l|2.51|l . The 
following relation is also useful: 

2?7s(k-p) = r) s (k)r] s (p) + r] d (k)r] d (p) 

+r) x (k) Vx (p) + r) y (k)i7„ (p) . ( A8) 



The facto r ?7 s (k- p) appearing in Eqs. I|A2|I . (|A3j) . (|A5jl 
and l)A6|l can be decoupled in this way. Even if only the 
part of r = s in Eqs. I|2.13[) and (|2.14|l is considered, 
there appear contributions from and p waves. 
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